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ABSTRACT 

Context. The dynamics of the solar core cannot be properly constrained through the analysis of acoustic oscillation modes. Gravity modes are 
necessary to understand the structure and dynamics of the deepest layers of the Sun. Through recent progresses on the observation of these 
modes - both individually and collectively - new information could be available to contribute to inferring the rotation profile down inside the 
nuclear burning core. 

Aims. To see the sensitivity of gravity modes to the rotation of the solar core. We analyze the influence of adding the splitting of one and several 
g modes to the data sets used in helioseismic numerical inversions. We look for constraints on the uncertainties required in the observations in 
order to improve the derived core rotation profile. 

Methods. We compute forward problems obtaining three artificial sets of splittings derived for three rotation profiles: a rigid profile taken 
as a reference, a step-like and a smoother profiles with higher rates in the core. We compute inversions based on Regularized Least-Squares 
methodology (RLS) for both artificial data with real error bars and real data. Several sets of data are used: first we invert only p modes, then 
we add one and several g modes to which different values of observational uncertainties (75 and 7.5 nHz) are attributed. For the real data, we 
include g-mode candidate, £=2, n=-3 with several splittings and associated uncertainties. 

Results. We show that the introduction of one g mode in artificial data improves the rate in the solar core and give an idea on the tendency 
of the rotation profile. The addition of more g modes gives more accuracy to the inversions and stabilize them. The inversion of real data 
with the g-mode candidate gives a rotation profile that remains unchanged down to 0.2 R Q , whatever value of splitting we attribute to the g mode. 
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1. Introduction 

The Sun is a magnetic star and it is now well recognized that the 
dynamical processes occurring in the solar interior are linked 
to the activity of the visible external layers (For a review of so- 



of the solar interior 



lar and stellar activity see Schrijver & Zwaan (2000)). The fact 
that the Sun is still active today, even at the present low rotation 
rate (in comparison to young stars), implies that the initial mag- 
netic fields are maintained, regenerated or amplified through 
dynamo effects which are induced by fluid motions within the 
star, namely rotation, convection and/or meridional circulation 



(support for this comes from observations, e.g., Hartmann & 



Noyes (T987JI). In this sense, it is necessary to reconstruct 



the solar internal rotational profile from the surface down to 
the core, to properly understand the magnetic activity of the 
Sun. Over the past decade, increasingly accurate helioseismic 
observations from ground-based and space-based instruments 
have given us a reasonably good description of the dynamics 
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Schou et al.|[T998l |Antia & Basul|2000l 



Thompson et al. 2003 and references therein). Helioseismic in- 
ferences have confirmed that the differential rotation observed 
at the surface persists throughout the convection zone. There 
appears to be very little, if any, variation of the rotation rate 
with latitude in the outer radiative zone (0.4 > r/R G > 0.7). 
The rotation rate is almost constant (« 430 nHz) in this region 
which is separated from the region of differential rotation by a 
narrow shear layer — known as the tachocline ( Spiegel & Zahn 
[T992l|Corbard et al.|1998] > — . 

The rotation profile of the Sun is also connected to different 
aspects of the structure and dynamics of the star. This is the rea- 
son why the rotation rate is needed to estimate the circulation 
and shear instabilities which are responsible for the redistribu- 
tion of chemical elements (Thompson et al.|2003| >. Moreover, 
the redistribution of angular momentum through the coupling 
between the turbulent convection and the rotation contributes 
to the strong differential rotation in the convective zone (e.g. 



Schou et al. 1998) and hence, to the dynamo effect that is 
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thought to be responsible for the 11 -year activity cycle ( |Brun| 
|et al.|2004{ |Gilman et al. 12007] and references therein) and the 
evolution of the Hale solar cycle ( Dikpati & Gilmanp 006 ). 

Although the helioseismic inferences in the radiative zone 
are not as precise as those found in the convective region, it 
could be confirmed that the rotation rate is flat and rigid down 
to approximately 0.4 R G . At least three processes have been 
proposed to explain this flatness of the rotation profile in the 
radiative zone but without great success: the redistribution of 
angular momentum by the effect of differential rotation which 



does not succeed in producing a completely flat profile ( Talon 



|& Zahn|1997| >; the effect of some magnetic fossil field instabil- 
ities that could flat the profile but they have not yet been found 
( Spruit 2002 ; Eggenberge r et al.|2005 1 and finally, some trans- 



port of angular momentum by internal gravity waves (Talon 



et al. 2002; Charbonnel & Talon 2005; Turck-Chieze & Talon 



2007 ). A general formalism has been developed recently to take 



into account all these different processes (Mathis & Zahn 2004 



2005), but more accurate observations of the solar rotation pro- 



file are needed to constrain the theoretical picture. 




Fig. 1. Rotational kernels for a p mode {( = 2, n=6) and the 
g-mode candidate (i = 2, n=-3). 



The analysis of the rotation profile at deeper layers (be- 
low 0.4 R Q ) and therefore, inside the solar burning core (where 
more than half of the solar mass is concentrated) could only be 
carried out with a few tens of p modes - the low-degree modes 
{( <3) -. Indeed, since the dawn of helioseismology when the 
works by |Claverie et al.| ( |198T]|1982| i led to the conclusion that 
the solar core rotates from 2 to 9 times faster than the surface 
rate, several groups have published different estimations of the 
rotation rate in the solar core - using acoustic modes - with 
contradictory results ( Jimenez et al.|1994||Elsworth et al.|1995| 
Fossat et al.|[l995l [Lazrek et al.||1996| |Chaplin et al.||2001| i' 
The importance of the low order p modes (below 2.3 mHz) 



to properly establish the profile below 0.4 R Q has been shown 



( Couvidat et al.|[2~003 i. Gravity modes having large sensitivi- 
ties to the solar core, will significantly contribute to establish 



the actual dynamical conditions of the core. This is illustrated 
in Fig.[T]representing the rotational kernel for a g mode (( = 2, 
n = -3) and for a p mode (( = 2, n — 6). This g mode is mostly 
sensitive to the region below 0.2 R Q whereas the p mode is sen- 
sitive to the region above 0.4 R G . It shows the importance of g 
modes compared to p modes for having access to the rotation 
of the core. 

The advent of the new millennium saw the burgeoning of 
the g-mode research based on the quality and accumulation 



of helioseismic data. In 2000, Appourchaux et aL] looked for 
individual spikes above 150 //Hz in the power spectrum with 
more than 90% confidence level that the signal was not pure 
noise. Although they could not identify any g-mode signature, 
an upper limit of their amplitudes could be established: at 200 
//Hz, they would fall below 10 mms in velocity, and below 



0.5 parts per million in intensity. Later, in 2002, |Gabriel et al. 



using the same statistical approach, found a peak that could be 
interpreted as one component of the £=l, n=\ mixed mode. 

A different approach based on the search of multiplets and 
recurrent signals in time (|Garcia & Turck-Chieze 1997 Palle 



& Garcia 1 1997[ |Turck-Chieze et al.||1998|l have been applied 
f[|JSoHC|^J velocity time series (Garcia et al. 



to GOLI 



2005). 



Some time-coherent patterns were found in the signal ( Gabriel | 
1999}, thus they could be potentially considered as g 



et al. 



modes. Turck-Chieze et al. (2004]) applied this technique to 
high-frequency multiplets in hope of reducing the detection 
threshold while maintaining the same confidence level. These 
authors found several patterns attributed to g-mode signals and, 
in particular, one was considered as a candidate for the mode 
t — 2, n — -3. In fact, |Cox & Guz ik (2004 ) postulated theoreti- 
cally that this mode could be the one with the largest amplitude 
at the solar surface. This candidate is still present in the anal- 
ysis of longer time series (Mathur et al. 2007, and references 
therein). Finally, the measurement of a signal that could be at- 
tributed to the separation in period of the dipole gravity modes 
and the comparison with solar models fosters a faster rotation 
rate in the core than the rest of the radiative zone ( jGarcfa et al.| 
[2007] ). 

In this work, we will study how the inferences about the so- 
lar core rotation profile could be improved by including gravity 
modes. We will study the effect of adding either one (the can- 
didate i-2 n=-3) or several g modes in the data set that will be 
inverted to infer the rotational profile. The effect of the obser- 
vational uncertainties on the derived rotational rate will also be 
analyzed, as well as the introduction of the g-mode candidate 
in real p-mode data sets. 



2. Methodology. 

Helioseismic inferences on the internal rotation rate of the Sun 
are carried out through numerical inversions of the functional 
form of the perturbation in frequency, Av„( m , induced by the ro- 



Global Oscillations at Low Frequenc y |Gabriel et al.|l995 
Solar and Heliospheric Observatory jDomingo et 



al.Jl995) 
al.fl995| ) 
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tation of the Sun, H(r, 6) and given by (see derivation in Hansen 
|et al.|1977> : 



Av„ f „, = - 1 - f f K ntm (r, ff)Q(r, 0)drd6 + e ntm 
2tt Jo Jo 



(1) 



The perturbation in frequency, Av„( m with error e„( m , that 
corresponds to the rotational component of the frequency split- 
tings, is given by the integral of the product of a sensitivity 
function, or kernel, K„( m (r, 6) with the rotation rate, Q.(r, 8), 
over the radius, r, and the co-latitude, 8. The kernels, K„( m (r, 8), 
are known functions of solar models. 

Equation [T] defines the forward problem for the solar in- 
terior rotation rate through global helioseismology, since it 
is possible to calculate estimates of the frequency splittings, 
Av„f m , that correspond to a given solar rotation rate, Q(r, 8). 



A is the matrix with the kernels, of dimension M x N and e is 
the vector containing the errors in D. 

The RLS solution is the one that minimizes the quadratic 
difference^- 2 = \Ax - D\ 2 , with a constraint given by a smooth- 
ing matrix, H, introduced in order to lift the singular nature of 



the problem (see, for instance, Eff-Darwich & Perez Hernandez 



1997). Hence, the function x is approximated by 



x est = (A T A + AH) 



l A T D 



(3) 



where A is a vector defining how much regularization is applied 
to each point x, of the inversion mesh, as introduced in |Eff-| 
|Darwich & Perez Hernandez | ( fl99"7] l. 

As a by-product of the inversion methodology, we could 
replace D from Eq. [2]to obtain 



(A T A + hHf 1 A T Ax = f Rx 



(4) 



Rigid 
Sfep 
Smooth 
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Radius (r/R ) 



hence 

R = (A T A + AHy l A T A 



(5) 



The matrix R, that combines forward and inverse mapping, is 
referred to as the resolution or sensitivity matrix, while the 
row is referred as the resolution kernel for the estimation of 
X; ( |Eff-Darwich et aL 2008). The diagonal elements R ti state 
how much of the information is saved in the model estimate 
and may be interpreted as the resolvability or sensitivity of x,. 
In this sense, it could be possible to use /?,-,• to see the effect of 
modifying the mode sets used in the inversion on the resolv- 
ability of each point of the inversion mesh. 



Fig. 2. Artificial rotation profiles for the solar interior as ex- 
plained in the text and used in the computation of the artificial 
data sets. The three artificial profiles have the same behavior in 
the convective zone. They incorporate latitudinal variations in 
the convection zone to mimic the actual rotation profile of the 
Sun. The plotted latitudes are 0° (solid line), 30° (dashed line), 
60° (dotted line) and 75" (dashed-dotted line). 
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The latter equation also defines a classical inverse problem 
for the sun's rotation. The inversion of this set of M integral 
equations - one for each measured Av„c m - allows us to in- 
fer the rotation rate profile as a function of radius and latitude 
from a set of observed rotational frequency splittings (hereafter 
referred to as splittings). The inversion method we have used 
is based on the regularized least-squares methodology (RLS). 
The RLS method requires the discretization of the integral re- 
lation to be inverted. In our case, Eq. 1 is transformed into a 
matrix relation 




Fig. 3. Observational error bars of p-mode splittings for degrees 
{=1,2, 10 and 20 as a function of the central frequency of the 
mode. Each degree is represented by a symbol as explained in 
the legend of the figure. 



D = Ax + e 



(2) 



where D is the data vector, with elements Av n [„, and dimension 
M, x is the solution vector to be determined at tabular points, 



A theoretical analysis was carried out in order to determine 
the effect of the addition of g modes on the derivations of the 
solar rotation rate of the burning core. Different artificial data 
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sets have been calculated using Eq. [T] and three artificial rota- 
tion rates D(r, 9) that are shown in Fig|2] They all have a dif- 
ferential rotation in the convection zone and a rigid rotation 
from 0.7 down to 0.2 R G equal to fi rz =433 nHz. In the first pro- 
file - the rigid profile -, that is our reference profile, the flat 
and rigid rotation includes the core. The second profile - the 
step profile -, is a step-like profile having a rate 3 times larger 
than the rest of the radiative zone below 0. 1 R Q and a rate of 
350 nHz in the region 0.1-0.2 R Q . Though this profile has un- 
realistic steep changes, it is useful to check the quality of the 
inversion as these steep profiles are difficult to reproduce. The 
rotational rate for the third profile - the smooth profile - in- 
creases gradually from 433 nHz at 0.2 R G reaching 1800 nHz 
in the centre, being in this sense compatible with the latest the- 
oretical studies. 



2000 



1500 



« 1000 



500 



0.4 




since it is more difficult to estimate the splittings as a conse- 
quence of the blending between the multiplet components of 
the modes due to the reduction in their life times (Bertello et al. 



2000 



et al. 



GarcfaetaLj200T] [2004{ |Chaplin et alpX)^ |Couvidat 



2003). Up to eight different g modes have been used 



in this work, four l=\ (with frequencies down to 100 //Hz) 
and four 1-2 (with frequencies down to 150 //Hz) which are 



the modes with the highest predicted amplitudes ( |Kumar et al. 
1996 Provost et al.||2000| ). Since g modes have not yet been 
characterized, different theoretical uncertainties have been used 
during the inversion process. Indeed the first value of 75 nHz 
for the uncertainty on the g-mode splitting, corresponds to the 



0.2 

Radius (r/R s ) 



tolerance in the search algorithm used by Turck-Chieze et al. 
(2004J and is related to a possible shift of the multiplet com- 
ponents due to the presence of a central magnetic field. Rashba 
et al.| ( |2007| l have already shown that this would shift the cen- 
tral frequencies of g modes in this region by such amount. The 
other value of 7.5 nHz is a typical uncertainty that could be ob- 
tained by fitting dipolar modes with ~ 4 years of data and could 
be a good example of what we could measure in the near future 
with the next generation of instruments. 

3. Results 

3.1. Inversions of artificial data 

A set of numerical inversions were carried out to study the ef- 
fect of adding g modes on the derivation of the rotation rate 
of the solar core. The analysis of the inversion results were 
complemented with the study of the resolution kernels of the 
inversions and the direct comparison of the sets of frequency 
splitting used. 



Fig. 4. Equatorial rotation profiles below 0.4 R e reconstructed 
with the p modes (Set Di) for the rigid profile (dotted line), the 
step profile (triple dotted-dashed line) and the smooth profile 
(dashed line). For the sake of clarity we have plotted the er- 
ror bars in the step profile only. The continuous blue, red and 
small dashed lines are respectively the step, smooth and rigid 
artificial rotation profiles. 



The different artificial data sets correspond to different 
mode sets, as explained in Table [T] The observational un- 
certainties for p modes (see Fig. (STT were calculated through 
Principal Component Analysis of the mode s ets extracted from 
a sample of 728 days-long MD^jtime-series ( |Korzennik|2"005l ) 
for p modes with degrees ranging from I -A to 25, whereas for 
t=\, 3 modes, the uncertainties were extracted from a com- 
bined GOLF-MDI time series ( |Garcia et al.|2004| >. The degree 
range of all data sets spans from l-\ to 25, however the fre- 
quency range of the artificial data sets depends on the degree 
of the mode, ranging from 1 to 2.3 mHz for (=1 to 3 and from 
1 to 3.9 mHz for (=4 to 25. As it is illustrated in Fig. [3] the un- 
certainties above 2.3 mHz for low degree modes are very large, 



j •! Olf «q «/J <>*, »4 »4: »a; 
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Al=3 



2000 2500 



Michelson Doppler Imager ( Scherrer et al. 1995 1 



Fig. 5. Difference of the p-mode splittings between the step 
profile and on the one hand the rigid profile (top) and on the 
other hand, the smooth profile (bottom). 



The inversion of the available p-mode splittings (see Fig|4]), 
as those included in set Di, reveals that it is not possible to re- 
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Table 1. Description of the artificial data sets used to study the sensitivity of p and g modes to the dynamics of the solar core. 



Freq. range (mHz) 



Data set 




g modes 






Uncertainty on g modes (nHz) 


p modes I = 1 , 3 


p modes I > 3 


Set Di 












1 < v < 2.3 


1 < y < 3.9 


Set D 2 


I 


= 2,n = -3 






75 


1 < v < 2.3 


1 < y < 3.9 


Set D 3 


f 


= 2,n = -3 






7.5 


1 < v <2.3 


1 < y < 3.9 


Set D 4 


£=l,n = -2 to 


-5 and I = 2, n = 


-3 to 


-6 


75 


1 < y < 2.3 


1 < y < 3.9 


Set D 5 


I = 1, n = -2 to 


-5 and ( = 2, n - 


-3 to 


-6 


7.5 


1 <y<2.3 


1 < y < 3.9 



cover any of the three artificial rotation profiles below 0.2 R Q 
(see Fig|2|. This result is also illustrated by the comparison of 
the splittings calculated from the rigid, step and smooth pro- 
files (see Fig|5]), since such differences fall below 1 nHz, be- 
ing the present level of uncertainties for these splittings above 
this value. The resolution kernels for these inversions (Fig|6]i 
also confirm the lack of sensitivity below 0.2 R Q , since it is not 
possible to properly locate and recover the resolution kernels 
below 0.2 Rq. 




0.4 0.6 
Radius (r/R Q ) 



the g mode (with an error bar of 75 nHz), when compared to 
the same resolution kernel calculated from the inversion of p 
modes. Unlike the splittings calculated for p modes, the differ- 
ences in the frequency splittings calculated from the three artifi- 
cial rotational profiles for the g modes could be larger than 200 



nHz (see Fig 10 1. It is shown in this figure that the differences 
of splittings calculated for a rigid profile and for the other two 
simulated rotation profiles, are around 200 nHz. This means 
that following the usual criteria of ~3 cr to have a proper detec- 
tion, a difference of 200 nHz is visible with the modes having 
an uncertainty of 75 nHz. Therefore, in this condition, it is pos- 
sible to discriminate between the rigid profile and the other one 
with a higher rotation rate in the core. 




Fig. 6. Resolution kernels computed in the inversion of the set 
D\ containing only p modes and calculated at two radii: 0.08 
Rq (solid line) and 0.16 R G (dashed line). 



When one g mode {1—2, n=-3 around 220 pHz) is added to 
the p-mode data set, as in the case of sets D2 and D3, the in- 
version results improve below 0.1 R Q (see Fig|8]and Fig J9]>, but 
there is not substantial improvement around 0.2 R Q . The match 
between the artificial rotational profiles and the profiles esti- 
mated from the inversions improves when the error assigned to 
the g mode is reduced. Unlike the case of the inversion of only p 
modes, the resolution kernels at 0.08 R G significantly improves 
when adding one g mode (see Pig j7j, in particular when the ob- 
servational uncertainty falls to 7.5 nHz. The resolution kernel 
at 0.16 R does not change (or slightly) with the addition of 



-1.0 I 1 1 1 I 1 1 1 I 1 1 1 I 1 1 1 I 1 1 1 L 

0.0 0.2 0.4 0.6 0.8 1.0 

Radius (r/R„) 



Fig. 7. Resolution kernels computed in the inversion of the set 
Di (i.e. including the g mode C=2 n=-3 with an error bar of 
75 nHz), D^(i.e., including the g mode with an error bar of 7.5 
nHz), D4 (i.e., including eight g modes with an error bar of 75 
nHz) and Ds (i.e., including eight g modes with an error bar of 
7.5 nHz) at 0.08 R @ . 



If eight g modes are added to the p-mode data set, as in 
sets D4 and D 5 , the inversion results (see Figs (8] and [9]) in both 
the level of uncertainties of the estimates and the matching to 
the proxy rotation profiles are significantly better than those 
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Fig. 8. Equatorial rotation profiles below 0.4 R Q reconstructed 
from the set Di (i.e., including the g mode 6=2 n=-3 with an 
error bar of 75 nHz), Z>}(i.e., including the g mode with an error 
bar of 7.5 nHz), D4 (i.e., including eight g modes with an error 
bar of 75 nHz) and D5 (i.e., including eight g modes with an 
error bar of 7.5 nHz) for the step profile. 
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Fig. 9. Equatorial rotation profiles below 0.4 R G reconstructed 
from the set D2 (i.e., including the g mode 6=2 n=-3 with an 
error bar of 75 nHz), Dj,(i.e., including the g mode with an error 
bar of 7.5 nHz), D4 (i.e., including eight g modes with an error 
bar of 75 nHz) and D5 (i.e., including eight g modes with an 
error bar of 7.5 nHz) for the smooth profile. Same legend as 
Fig§ 



obtained from the inversion with only one g mode. The dif- 
ferences in the frequency splittings calculated from the three 
artificial rotational profiles for all the eight g modes are signif- 
icantly larger than the observational uncertainties and hence, 
new information could be gained as compared to that given by 
just one g mode. This is particularly important in the presence 



of noisy data, since the larger the number of g modes, the better 
the averaging of the unwanted effects of the noise in the data 
will be. The addition of several g modes helps to better define 
the resolution kernels below 0. 1 R (see the resolution kernel 
at 0.08 R Q Fig. [7j, whereas the resolution kernel at 0.16 R G 
does not significantly change by the addition of g modes. In this 
sense, the inferences about the rotational rate of the core will 
be significantly improved below 0.1 R Q , where the energy of 
the g modes is maximum. Very high frequency p modes (above 
2.5 mHz) for 6=1 and 2 should be characterized to better define 
the region between 0.15 and 0.25 R Q ( Garcia et al.|2008) 




Fig. 10. Splittings differences of g modes between the rigid pro- 
file and on one hand, the step profile and on the other hand, 
the smooth profile (top). We have only represented the eight 
g modes used in the inversions. The filled square corresponds 
to the g-mode candidate (=2 n=-3. We have also plotted the 
splittings of these modes for the rigid profile (bottom). 



3.2. Inversion of real data 

We have used inversions to study the compatibility of present 
p-mode frequency splittings with the splittings estimated by 
Turck-Chieze et al. (2004) for the £=2, n=-3 g mode. The p- 
mode set corresponds to the mode set used in D 1; where the 



splittings correspond to those calculated by Korzennik (2005) 
for 2088 day-long MDI time-series. In Turck-Chieze et al. 



(2004), three scenarios were proposed to explain the detected 
pattern around 220 //Hz, with two possible values for the split- 
tings, namely 300 nHz, if this were a detection of two modes (a 
combination of an 1=2 and a ^=5), and 600 nHz, if all the vis- 
ible components correspond to the same mode (the 1=2 n=-3 
which implies an inclined core rotation axis). 

Five inversions are carried out (see FigfTT}, namely one in- 
version including only p modes and other four containing one 
g mode, but with different estimates of the frequency splitting 
(300 and 600 nHz) and two observational uncertainties (75 and 
7.5 nHz). As it was illustrated in Figs. |4j [5] and [6} there is no 
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sensitivity of the observed p modes to the dynamics of the in- 4. Conclusions 



ner solar core. Different rotational profiles below 0.2 R Q are 
obtained for the different combinations of the value of the in- 
troduced g-mode splitting and its corresponding uncertainty. In 
all cases, these values are compatible with the data calculated 
for the p modes in the sense that the inversions are unchanged 
and stable (e.g. the inversion does not show any oscillatory be- 
haviour) above 0.2 R Q when the g mode is added to the data 
set. In Fig j8] an oscillation around 0.2 R Q appears with the ar- 
tificial data when g modes are associated to an error bar of 7.5 
nHz. This is also observed with the real data. This is an artefact 
of the inversion. The rotation profile obtained using the high- 
est value of the g-mode splitting (600 +7.5 nHz) proposed by 



Turck-Chieze et al. ( 2004 1, gives a rate in the inner core that is 
compatible with the result from the dipole analysis carried out 
by |GarciaetaT1 ( p007l l. 



To quantify the compatibility between the p-mode data and 
the g-mode candidate, we have calculated the normalised resid- 
uals for all the observed modes (C,m,n) defined by (errata - 
<5vinv)/o% where 5vd a ta is the value of the splitting in the data, 
6v mw , the value corresponding to the rotation profile obtained 
with the inversion and cr, the error bar associated to the split- 
ting of the mode. Table [2] gives the mean value of these resid- 
uals for the low-degree p modes {( <3, below 2 mHz) and the 
g-mode candidate for the four inversions (the two values for 
the splittings and the two values for the uncertainty for the g 
mode). We can see that the difference of splitting between the 
real data and the results of the forward problem on the inferred 
rotation profile is less than 1.5 cr for the p modes. However, 
for the g-mode candidate, this difference goes up to ~3 cr. This 
is due to the fact that the rotation profile has some uncertain- 
ties that have an impact on the splittings calculated. Globally, 
the results with the g-mode candidate are compatible with the 
information contained in the observed p-mode data. 



In this paper, we have studied how the inversion of several ar- 
tificial rotation profiles could be improved when g modes are 
added to the present set of observed p modes. The introduc- 
tion of one g mode -the candidate 1-2, n=-3- significantly im- 
proves the solution in the inner core (below 0. 1 R & ), when com- 
pared to that obtained using only p modes. It gives the general 
trend of the solar core rotation but there is not accurate infor- 
mation on the profile itself. If more g modes are added to the 
inversion data set (four t-\ and four sectoral <"=2), the result in 
terms of accuracy and error propagation, improves compared 
with the inversion including only 1 g mode. However there is 
still information missing in the region between 0.1 - 0.2 R Q , 
where the energy of the g modes is significantly lower than in 
the region below 0.1 R . The information given by the p modes 
is negligible due to the lack of sensitivity to these depths, the 
high level of uncertainties we have in their determination and 
the noise present in the data. 

Finally, for the real data, the rotation profile obtained us- 
ing the highest value of the g-mode splitting gives a rate in 
the inner core that is compatible with the result obtained with 
an independent technique by Garcia et al. ( 2007| l, if we put an 
error bar of 7.5 nHz. Moreover, we obtained a limit down to 
which we can trust the inversion of the real data. All the values 
proposed for the splittings of the £-2, n=-3 g-mode candidate 
are compatible with the splittings calculated for the p modes. 
Indeed, having in mind that the small oscillation is related to 
the inversion and not to the data, the addition of the g-mode 
candidate with different values for the splitting and their uncer- 
tainty, does not change the estimated profile above 0.2 R Q . 
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Fig. 11. Equatorial rotation profiles reconstructed with the real 
data as explained in the text and adding the g-mode candidate 
1=2 n=-3 with two different splittings (300 and 600 nHz) com- 
bined with two different error bars (75 and 7.5 nHz). 
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